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To explain the observed decay of superfluid turbulence at very low temperature, it has been 
proposed that a cascade of Kelvin waves (analogous to the classical Kolmogorov cascade) transfers 
kinetic energy to length scales which are small enough that sound can be radiated away. We report 
results of numerical simulations of the interaction of quantized vortex filaments. We observe the 
development of the Kelvin-waves cascade, and compute the statistics of the curvature, the amplitude 
spectrum (which we compare with competing theories) and the fractal dimension. 

o 

CN ; PACS numbers: 67.25.dk, 47.37.+q 03.75.Kk 

I. MOTIVATION 

■ ^ ■ Quantum turbulence [H-Q consists of a disordered tangle of reconnecting superfluid vortex filaments. Because of 
O ' quantum mechanical constraints on the rotational motion, these filaments carry the same quantized (fixed) circulation 
r and are very thin: the (fixed) vortex core radius, ao, is much smaller than the average distance between the filaments, 
I ' £. Quantum turbulence is easily created by agitating superfluid helium ("^He) with propellers grids 0, forks @ 
"j^ , or wires 0, by applying a heat flow [8l-[l2j or by injecting a stream of ions It is also studied in superfluid ^He-B 
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[3, [3] and, more recently, in atomic Bose-Einstein condensates [3, E 
If the temperature is relatively large (more than 1 K in ^He) , the turbulent kinetic energy contained in the superfluid 
d • vortices is transferred by the mutual friction [isj into the viscous gas of thermal excitations (the normal fluid) and then 
decays into heat; therefore a constant supply of energy (continuous stirring for example) is needed to maintain the 
intensity of the turbulence. If the temperature is relatively small (less than 1 K in '*He), the normal fluid is ne gligi ble, 
■ O ' but, despite the absence of viscous dissipation, the turbulence still decays [l^, [lO]- The Kelvin- waves cascade j2ll - [23l | 
^ , was proposed to explain this surprising effect. 

A Kelvin wave is a rotating sinusoidal or helical displacement of the core of vortex fllaments away from its unper- 
turbed position l24[-[2q . The dispersion relation of a Kelvin wave of angular frequency w and wavenumber k along a 
straight vortex is [27| 



> 

\^ ■ where ao is the vortex core radius and Kn is the modifled Bessel function of order n. In the long-wavelength 
' approximation (fcap ^1), the angular frequency reduces to 

o ■ 

TP TP 
~ (ln2/(fcao) - 7) = (lnl/(fcao) - 0.116) (2) 

•T— I , 

■ where 7 = 0.5772 is Euler's constant; the negative sign in this expression means that the Kelvin wave propagates in 
^—^ ' the direction opposite to the orientation of the unperturbed vorticity. Eq. [T]was originally derived for a thin, hollow- 
cored vortex in a perfect Euler fluid; a similar expression was obtained for a vortex in a Bose-Einstein condensate 
[2^ [2^. The simplified dispersion formula uj{k) = cfc^, where c is a constant, is often used in the Kelvin- waves cascade 
literature. 

The Kelvin- waves cascade is the process in which the nonlinear interaction of Kelvin waves creates waves of shorter 
and shorter wavelength A = 2iT/k. At high enough temperature the mutual friction would quickly damp out the 
shorter Kelvin waves [l^l, but at low temperatures the cascade proceeds unhindered, until k is large enough that 
sound is efficiently radiated away (phonon emission) by rapidly rotating vortices [3(]l - l33l |. There is thus an analogy 
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between the Kelvin-waves cascade in quantum turbulence (where the energy sink is acoustic) and the Kolmogorov 
cascade [H, in classical turbulence (where the energy sink is viscous) . 

It is thought that in superfluid helium at low temperatures both the Kolmogorov cascade and the Kelvin-waves 
cascades operate at the same time [s^. Let A = L/V be the vortex line density (vortex length L in the volume 
V) of a homogeneous isotropic vortex tangle, and £ « A~^/^ be the average intervortex spacing. The Kolmogorov 
cascade (which in a classical fluid is often described as the break up of eddies (37l | ) involves rearrangements or partial 
polarisationjs^ of vortex lines over length scales of wavenumbers k « 2tt / 1. The Kelvin- waves cascade occurs on 
individual vortex filaments for wavenumbers k » It: j I. Recent work has raised the possibility of the existence of a 
bottleneck [s^ between the two cascades, for which energy would pile up at some length scale of the order of I. 

Direct evidence of the Kolmogorov cascade in superfluid helium was provided by the experiment of Maurer and 
Tabeling [4l[. Numerical evidence is also convincing: the Kolmogorov ~ fc~^/'^ energy spectrum has been seen in 
calculations performed using both the Schwarz filament model [43 and the Gross-Pitaevskii equation j43l - l45j which 
describes a Bose-Einstein condensate. On the contrary, the Kelvin-waves cascade is less understood; there is no direct 
observation of it, and the theory is mired in a controversy between approach es prop osed by L'vov, Nazarenko and 
collaborators [s^, HO, EB-IEll on one hand, and by Kozik and Svistunov [111, [s^, I52l - l56l | on the other. The two theories 
predict similar Kelvin-waves amplitude spectra, but differ on matters of principles and in some important details. 

Lacking direct experimental evidence of the Kelvin-waves cascade and its properties, experimentalists currently 
attempt to identify indirect but observable effects of the cascade and of any bottleneck; for example they study the 
temperature dependence of the observed effective kinematic viscosity jsvl ISq. The results are still inconclusive. 

The aim of our work at this stage is not to review or compare rival theories, but to add information. We present 
results of numerical calculations of the Kelvin-waves cascade. The few published calculations differ in details (e.g. 
whether the initial vortex is curved or straight, whether the cascade is forced continuously or results from an initial 
condition - either an artificially imposed cusp or a dynamically resolved vortex reconnection) . More information is 
therefore necessary to learn the main physical features of the turbulent Kelvin-waves cascade and to assess which 
results are robust. In particular we shall focus the attention to the magnitude and spectrum of the Kelvin waves, the 
statistics of the vortex curvature and the fractal dimension. 

Finally we note that, in presenting our results, we shall pay particular attention to how well energy is conserved in 
a given numerical approximations. In the absence of mutual friction, vortex evolution in between vortex reconnection 
events and without a turbulent cascade which radiates sound, energy should be conserved; yet previously published 
numerical calculations of the evolution of superfluid vortex filaments do not report how well the conservation law is 
satisfied. 



II. MODEL 



Hereafter we use parameters which refer to superfluid ^He: circulation F = 9.97 x 10^"' cm^/s and vortex core 
radius ao ~ I0~^ cm, but our results can be generalised to turbulence in low temperature "^He-B. Following Schwarz 
[sgj . we describe vortex filaments as space curves s = s(^,t) where ^ is arc length and t is time. In the absence of 
mutual friction and of any externally applied superflow, the self-induced velocity of a superfluid filament at the point 
s is given by the Biot-Savart law 



ds F / (s-r) 

The line integral extends over the entire vortex configuration which is discretised into a large number of points 
(i = I, • • • N'). The singularity at s = r is removed in a standard way by considering local and non-local contributions 
to the integral. If is the position of the z*^ discretization point along the vortex line, Eq. [3] becomes [gO] 



Here £i and li+i are the arclengths of the curve between points Si_i and and between and s^+i, and C is the 
original vortex line without the section between Si_i and Si+i. For simplicity we assume that the cutoff parameter is 
a « ag. Other de-singularization techniques are also used in the fiuid dynamics literature [gH l62j. 

The number of discretization points, iV, changes with time. As the simulation progresses, new discretization points 
are introduced to maintain the resolution along the vortex filament. If the separation between two points, and 
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Sj_|_i, becomes greater than some threshold 5 (which we call the minimum resolution), we introduce a new point at 
position Sir given by 



- ^(s. + s.+i) + i^^m lil, Rej (5) 

where i?;' = |Sj, |~^. In this way s^, = (Sj +Sj^]^)/2, that is to say the insertion of new points preserves the curvature. 
Points are removed if their separation is less than S/2, to ensure that our shortest length-scale is fixed. The above 
algorithms ensure that the separation of vortex points along a filament lies between 6/2 and 6. 

For accuracy, we approximate all spatial derivatives using 4*^ order finite difference schemes which account for 
varying mesh sizes along the vortex filaments (63| . Let be the i^^ point on the vortex filament; the two points behind 
have positions Si_2 and Si_i, and the two points in front have positions s^+i and Sj_|_2- We denote = \si_i — Si„2|; 
ii = \s.i " Si_i|, li+i = \si+i - Sj|, and 1^+2 = |sj+2 - Sj+i|- We can then write 

s- = AiiSi^2 + BiiSi^i + CuSi + DiiSi+i + EuS,+2, (6) 

where the coefficients An, Bu, Cu, Du and En are given by, 

, ^ i^£l+,+e^i^+l£^+2 

^,_i(^,_i+£,)(^,-i+^. + ^.+i)(^.-i+^, + ^,+1+4+2) ^' 
Cu = -iAu + Bu + Du + Eu) (9) 

£,+2(^.+ l+^*+2)(^*+^ 

In a similar manner we can write, 

s" = A2,;Si_2 + i32iSi_i + C2iSi + D2iS.i+i + E2iSi+2, (12) 

where the coefficients A2i, i?2o C2ij D2i and i?2i arc given by. 



i+2\ 



4-1 (4-1 + 4)(4-i + + + 4 + 4+1 + 4+2) 



(13) 



^ 2[24-l4+l + 244+1 - iUl + ^^-l^^+2 + 44+2 - 4+l4+2] 

4-i4(4 + 4+i)(4 + 4+i+4+2) ^ ' 

C2^ = -[A2^ + ^2, + D2, + ^2,) (15) 

^ _ 2[-4-i4 - ^? + 4-i4+i + 244+1 + 4-i4+2 + 244+2] , . 

4+i4+2(4 + 4+i)(4-i+4 + 4+i) ^ ' 

^ ^ 2[4-i4 + ^f-4-i4+i- 244+1] , . 

" 4+2(4+1 +4+2)(4 + 4+1 +4+2)(4-i +4 + 4+1 +4+2) ^ ' 
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Note that if we set ii-i = £i = ii+i = £^+2 = h, then the above expressions reduce to famihar finite-difference schemes 
for a uniform mesh: 

= Y^(s»-2 - 8s,_i + 8s,+i - s,+2) + 0{h^), (18) 



s" = Y^(-s»-2 + 16s,_i - 30s, + 16s,+i - s,+2) + (19) 

The Biot-Savart law describes an incompressible flow. In order to model energy dissipation at large k (phonon 
emission), we also remove points if the local wavelength is smaller than a specified value. Since the smallest length- 
scale is (5/2, the maximum curvature is of the order of 2/(5. At each time step we compute the local curvature 
C(^) = |s"| = \cfs/d^'^\ at each discretization point Sj [j = 1, • • • N) If, at some location along a vortex filament, the 
local curvature exceeds the critical value 1.9/(5 cm^^ (95% of the maximum value 2/(5), points are removed and the 
filament is smoothed locally. This leads to a loss of line length, hence a loss of energy. 

It must be stressed that we are not interested in the phonon emission per se, but rather in the consequence of the 
presence of an energy sink on the Kelvin waves cascade. The actual numerical value of 5 is chosen for numerical 
convenience and is orders of magnitude larger than the microscopic lengthscales which are physically relevant to 
phonon emission and vortex reconnections. Therefore our model contains (and controls) a mechanism of dissipation, 
and is different in scope from vortex methods dedicated to perfect Euler dynamics. In the next section we describe 
results obtained with two different spatial resolutions: 5 = 10~^ cm (high resolution) and (5 = 5 x 10~^ cm (low 
resolution). We shall see that our results are independent of the resolution used. 

Time evolution is based on a order Adams-Bashforth scheme. Consider an evolution equation of the form 
ds/dt = V. We time-step the equation according to 

sr+' = ^7 + ^(23vr - i6vri + svr^) + o{At') (20) 

where At is the time step and the superscript n refers to the time t„ = nAt (n ~ 0, 1, 2, • • •). Lower order schemes 
are used for the initialisation, when older velocity values are not available. The time step At is chosen so that the 
period Tmin = ^Tr/ujrnax of the fastest Kelvin wave (which has wavenumber kmax — 27r/Amm ~ t^/S) contains at least 
25 time steps; we have verified that reducing our time step further gives the same results. 

If two discretization points become closer to each other than the local discretization distance, our numerical algo- 



rithm reconnects the two filaments [6J| subject to the criteria that the total length (as proxy for energy) decreases 
. In this way the same quantity 6 controls both energy dissipation mechanisms in our model - sound emission and 
reconnections. Self-reconnections (which can arise if a vortex filament has twisted and coiled by a large amount) are 
treated in the same way. Since reconnections involve only anti-parallel filaments, prior to reconnection we form local 
(unit) tangent vectors s' = ds/d^, where ^ is arc length, and, using the inner product, we check that the two filaments 
arc not parallel. Finally, We remove any small loops with less than three discretization points. 

During the evolution we monitor the length L and the kinetic energy E (in unit of the density) of the vortex 
configuration; the latter is 



E^\j^^,^dV. (21) 

If the calculation is performed in a periodic box, the vortices can extend across the boundaries; the only way to 
compute the energy integral is to discretise the integrand on a three-dimensional spatial grid, computing the velocity 
at each grid point; unfortunately this procedure converges too slowly to be practical. However, if the calculation is 
performed in an infinite volume, the vortex filaments are closed loops. Assuming that the velocity field decays to zero 
at infinity, the energy integral can be transformed into a line integral (66| 



E=-[ v^dV = / V • r X udV = r / V • r X s'd^, (22) 
2 Jv Jv J c 

which converges much better. In writing Eq. [22] we have used the fact that the vorticity a; = V x v is concentrated 
on the vortex filaments [ujdV = Ts'd^). 

In the absence of mutual friction, during Biot-Savart evolution (but not during vortex reconnections, and not while 
the Kelvin- waves cascade is in progress with the associated sound emission which we model), the kinetic energy E 
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is conserved, hence changes in E measure the accuracy of the numerical method. Our current method conserves 
energy about a factor of ten better than our preHminary unpublished calculations. The energy-conserving properties 
of previous calculations in the supcrfluid literature is unknown. 

Unlike the energy E, the vortex length L is not conserved. However, if the Biot-Savart law is replaced by the Local 
Induction Approximation (LIA), the length is conserved. The LIA was used in the original work of Schwarz [59| and 
by others to reduce the computational cost. We do not use the LIA in our calculations. 

A quantity which is particularly interesting is the normalized histogram (probability density function, of PDF for 
short) of the curvature C(^), sampled on the discretization points: we call this quantity PDF(C). 

III. RESULTS 
A. Cascading and non-cascading vortex configurations 

In the first set of numerical experiments we consider a single vortex filament set across a periodic box as in Fig. [1] 
The initial shape of the vortex is not straight, but is perturbed by two sinusoidal Kelvin waves of amplitude A and 
different wavelengths A. We compute the time evolution and observe that that this initial condition does not yield 
a cascade to shorter waves. However, if we start with three parallel vortex filaments (each filament with two Kelvin 
waves), as shown in Fig. [^Jlcft), we notice that, as the vortices rotate around each other, they become more wiggly, 
see Fig. [21[right). 

This result suggests that, besides vortex rcconncction events (as shown in Rcf. [l^), the Kelvin- waves cascade can 
be initiated, more in general, by the interaction of nearby curved vortices via the velocity field which they induce 
on each other. The result is also consistent with wave-turbulence theory [67l |. which states that, for one-dimensional 
waves obeying a dispersion law of the form uj{k) = ck^, a four- waves process is not sufficient to generate a cascade. 
According to the theory, the minimal initial condition for the cascade consists of three waves which satisfy the resonant 
condition that fci + fc2 + fcs = fc4 + fc5 + fcg and a;(fci) -1-^(^2) -|-a;(fc3) = ^(^4) -I- ^(fcs) -I- ^(fcg), where indices 1, 2, 3 and 
4,5,6 refer to incoming and outgoing waves respectively. The initial condition shown in Fig. [^l^left) contains more 
waves than the minimal number of waves (if it contained only the minimal number of waves we would need to seek 
the exact resonant condition, which is not our purpose at this stage). 

In order to pursue this investigation more quantitatively, we want to monitor the kinetic energy of the vortex 
configuration. Unfortunately, as we have mentioned in Section HIl if the vortex configurations is set in a periodic box, 
it is not possible to accurately determine its energy. Hereafter, in order to use Eq. I22[ we consider only closed vortex 
loops contained in an infinite volume. 

In a second set of numerical experiments we start with a single vortex ring in an infinite volume; the ring is seeded 
with two Kelvin waves, as shown in Fig. |31 During the time evolution we observe that the waves interact, but there 
is no cascade, consistently with 

the evolution of the vortex line of Fig. [1] Indeed, if we examine the PDF of the curvature, we notice that PDF(C) 
does not evolve in time to larger values of C, as shown in Fig. St^top). Unlike the vortex line in the periodic box of 
Fig. [TJ now we can monitor how well the kinetic energy is conserved in our numerical calculations. Fig. |3Kbottom) 
shows that the relative change of energy during the time evolution of this non-cascading vortex configuration of 
interacting Kelvin waves is only AE/E « ±0.06 percent. This value describes how well we can integrate Euler 
dynamics. 

B. Curvature and energy 

Proceeding in analogy with Fig. [51 we now consider a third set of numerical experiments with three vortex rings 
oriented in the same direction as shown in the first snapshot of Fig. [S] each ring is initially perturbed by two random 
Kelvin waves. This set-up contains enough waves to trigger the Kelvin-waves cascade and allows us to carefully 
monitor the energy. 

The initial total number of discretization points is A^'o — 600; at the end of the calculation at time t = 2 s we have 
N = 773. The minimum distance between discretization points is set to (5 = 5 x 10~^ cm. All three rings are oriented 
in the same way, hence they travel in the same (negative z) direction, leapfrogging around each other. This simple 
vortex bundle configuration evolves into a Kelvin- waves cascade: Fig.[S]shows that during the evolution the wiggliness 
of the vortices increases. This is in agreement with the straight vortices of Fig. [21 Since the wiggly vortices of Fig. [5l 
are relatively close to each other, there are many reconnections up to time t = 0.1 s, less reconnections between 
0.1 < < < 0.5 s, and no reconnection afterwards. The reconnections do not appear to affect the vortex configuration 
in a significant way: the vortex bundle seems structurally robust. 
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We start our analysis of the Kelvin-waves cascade by noticing that it induces a stretching of the vortices: Fig. [6] 
shows that the relative increase of vortex length is AL/L « 6%. This effect (the transformation of part of the 
interaction energy into length) is in agreement with previous calculations of reconnecting vortex bundles performed 
using the Gross-Pitaevskii model [sl, HI] . 

The behaviour of the energy E is shown in Fig. [T] During the initial stage of the evolution (0 < t < 0.5 s) 
the energy slightly declines, then it becomes approximately constant, with relative fluctuations AE/E of less than 
3%. The energy fluctuations are larger than for non-cascading vortex configurations, because in the cascading case 
the presence of short, rapid Kelvin waves induces rapid velocity fluctuations. Fig. Hljtop) and Fig. [SJbottom) show 
respectively the maximum speed Vmax and the maximum acceleration [dv / dt)max sampled over the discretization 
points Sj (j = 1, • • - iV). It is apparent that velocity fluctuations are larger during the initial stage of the cascade 
(0 < t < 0.5 s). We shall return to why the energy initially decreases after discussing the vortex curvature. 

Fig. ini shows the time evolution of mean, maximum and minimum curvature, < C >, Cmax and Cmim sampled 
over the discretization points Sj. Fig. [9jtop) shows that < C > undergoes a rapid initial rise, then slowly decays, 
almost settling down. Fig. [9][bottom) shows that the maximum curvature is always much less than the maximum 
value of curvature allowed by the numerical resolution, which is 1900 cm~^ and is represented by the dashed line. 
The curvature PDF is shown in Fig. (TU] Note that, as time progresses, PDF(C) moves to the right towards larger 
values of C. This effect must be compared against the much smaller movement to the right of the curvature PDF for 
the non-cascading vortex configuration which is showed in Fig. IH^top). Numerical experiments suggest that PDF(C) 
is indeed a convenient quantity to monitor if one looks for evidence of the Kelvin-waves cascade. At later stages, as 
energy is lost at small scales, we observe that the turbulence decays, as in the experiments [l^ll^l, and the curvature 
PDF moves to smaller values of C. 

It is also interesting to note in Fig. [TU] that, after an initial transient, the curvature PDF assumes a power law 
dependence PDF(C) ^ C" for large values of C. Fig. [TT] shows that the exponent a quickly settles down to the 
average value a ~ —3.5. In the low resolution run {S = 10"'^ cm) a similar calculation gives a = —3.34, which 
confirms that our results are independent of the numerical resolution. 

We can now return to the question as to why there is a small energy loss during the initial stage of the evolution 
from t ~ to t ~ 0.5 s in Fig. [T] Our numerical algorithm has two mechanisms through which length, hence energy, 
can be lost: vortex reconnections and smoothing of regions of high curvature. Which is more important ? To find 
the answer we repeat the calculation with the vortex reconnection algorithm switched off. The results arc shown in 
Fig. 1121 If we compare the energy computed with reconnections (Fig. [7]) and without reconnections (Fig. ll2f top)) we 
see that E is more constant without reconnections. It is also instructive to compare curvature statistics computed 
with reconnections (Fig. ^ and without reconnections (Fig. [121 bottom three graphs). It is apparent that vortex 
reconnections are the most important energy sink at large k: without reconnections, the maximum curvature Cmax 
is much larger, and often reaches the numerical limit represented by the horizontal dashed line in Fig. II 2f bottom). 
Since Cmm is approximately the same, < C > is also larger than in the reconnecting case. Therefore, without vortex 
reconnections the cascade proceeds to the smallest scale which is numerically allowed, at which point the smoothing 
algorithm kicks in. We conclude that the initial energy loss shown in Fig.[7|is due to the frequent reconnections which 
occur during the initial stage of the evolution. 



The next quantity of interest is the amplitude a(^) of the Kelvin waves and the spectrum of the Kelvin- waves 
cascade. Following Svistunov [2l| . we introduce the concept of smoothed vortex filament: at every time t, we use every 
n = 15 points as nodes of a cubic-spline interpolation to obtain a new (smoothed) filament, Sgmooth from the original 
discretised vortex filament Sj (j = 1, • • • N). Fig. [T3] (top) is an example of this procedure: in the figure we see the 
actual discretization points Sj {j ^ 1, • • ■ N) together with the smoothed line. We choose n = 15 because, in a range 
which is approximately 10 < n < 20, the amplitude spectrum (which we shall define shortly) docs not depend on n, 
as shown in Fig. I13f bottom), confirming that our numerical smoothing procedure is robust enough. 

We then define the Kelvin wave amplitude a(^) as the distance between the original filament and the smoothed 
filament: a(^) = |s — Sgmoothl- The mean amplitude < a > as a function of time is shown in Fig. 1141 Note that the 
mean amplitude < a > saturates for t > 1 s. 

Finally, the amplitude spectrum, A(k), is defined as 



and is shown in Fig. [151 We find that the amplitude spectrum behaves as A{k) ^ k^ for large k. The best fit is 



C. Spectrum 




(23) 
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/? = —3.10; we obtain /3 = —3.21 in the low resolution run. Fig. [TBI shows that the scaling exponent (3 is independent 
of time, as the Kelvin-waves cascade develops. 



D. Angle 

From the dot product of the local tangent vectors, s' and s^j^j^^^j^, at each time t we calculate the angle 6 between 
the original vortex filament s and the smoothed filament Ss,„ooth- The normalised histogram, PDF(0), computed at 
t ~ 0.002 s is shown in Fig. [TTf top). The peak value is approximately at w 15°. 

In order to compare the distribution which we find with that of a standard normal distribution (zero mean and 
unit standard deviation), we make use of the following standard result from statistics. All normal random variables 
with mean /i and variance a can be 'standardised' to fit a standard normal distribution by subtracting the mean and 
scaling this quantity by the standard deviation. Fig. ([T7]) (bottom) shows a plot of zq = {9 — 9)/a{9) where 9 and 
<7{9) are respectively the mean and the variance of the distribution of the angle 9. It is apparent that the distribution 
of the angles is close to a normal distribution. This result validates the weak-turbulence theory approach to the 
Kelvin- waves cascade problemjT^. A small departure from Gaussianity is present at large angles; it is likely that this 
effect is due to the fractal nature of the vortex filament (see section IIlIF[) . 



E. Reconnecting vortex rings 

The last set of numerical experiments is concerned with two vortex rings which reconnect with each other. This 
configuration is similar to the collision of four vortex rings which was used in a previous workj2^. Fig. 1181 shows the 
time evolution. It is apparent that, over the time scale which is considered, no Kelvin-waves cascade develops. This 
visual conclusion is confirmed by the inspection of PDF(C), shown in Fig. ll9f top). which does not progress in time to 
larger values of C, unlike Fig.[TOl Fig. [TOlbottom) shows that the energy remains constant within AE/E ^ ±0.12%. 

If the initial vortex rings are not perfectly circular, however, the evolution proceeds toward a Kelvin-waves cascade, 
as shown in Fig. [201 The increasing wiggliness of Fig. [501 must be compared against Fig. [THl Following the initial 
reconnection of the two rings, we observe 4 more reconnections, which are actually self-reconnections since there is only 
one vortex filament left. Fig. I20f c) shows a self-reconnection event which arises from a large twisting (supercoiling) of 
the vortex filament, and which results in the emission of a small vortex loop. This self-crossing scenario was predicted 
by Kozik and Svistunov[55|. Similar small- loop emission were also observed in numerical studies of vortex bundle 
dynamics [68| . 

The PDF of the curvature confirms the existence of a cascade: Fig. [HJtop) shows that PDF(C) progresses in time 
to larger values of C. Fig. [?IT bottom) shows the behaviour of the energy, which has the same slight initial decreases 
as in Fig. [3 which we have discussed. 



F. Fractal dimension 



To quantify the wiggliness of the vortex configuration undergoing the Kelvin-waves cascade, we compute the correla- 



tion dimension introduced by Grassberger and Procaccia |7l| . The correlation dimension provides an upper-bound 
estimate of the fractal dimension The definition is the following: given the N discretization points (j = 1, • • • N), 
we define the correlation integral as K(e) = n/N"^ where n is the number of points whose separation is less than e. In 
the limit of e 0, K{t) takes the form, 

K(e)^e^-, (24) 

where cd is the correlation dimension. This quantity can be calculated efficiently at the same time as the Biot-Savart 
integral is calculated. The result is shown in Fig. [52] for reconnecting vortex rings undergoing the Kelvin-waves 
cascade. Note the rapid increase of cd during the initial stage caused by reconnection events, followed by a decline 
due to the numerical dissipation at high wavenumbers which models physical dissipation (sound emission) . 



IV. DISCUSSION 



Our numerical calculations show that the Kelvin- waves cascade can be initiated by the interaction of neighbouring 
vortex lines as well as by vortex reconnections. We find that, shortly after an initial transient, the Kelvin- waves 
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cascade quickly settles in and saturates: Fig. (Illf) and Fig. p6)) show that the slopes a and /? of the curvature 
histogram at large C, PDF(C) ~ C", and of the amplitude spectrum at large fc, A{k) ~ fc^, remain approximately 
constant. 

We also find that the Kelvin-waves amplitude spectrum scales approximately as A(fe) ~ k~^'^ for large fc, in 
agreement with A ^ reported in Ref. [2^ for a continuousl y e xcited vortex line, and less steep but in fair 
agreement with the theories of L'vov, Nazarenko and collaborators jSlj and of Kozik and Svistunov [sB], who predict 
n-k ~ = and ^ fc^^^/^ — k~^ '^ respectively, where their kelvon occupation number Uk is equivalent to 

our A{k). The applicability of wave turbulence theory to the problem is supported by our finding that the angle 9 
between the smoothed vortex line and the actual line is small and has approximately a Gaussian nature, see Fig. [171 

The regime which we have identified in our calculation is clearly a weak turbulent regime, as demonstrated by 
the small angles between the smoothed vortex filament and the actual filament - see Fig. fTWtop). We think that 
in a stronger regime larger twists would appear and induce more self-reconnections than the few which we observe 
in Fig. [20l In must be stressed that our numerical cutoff effectively short-circuits the formation of loops at scales 
smaller than 6. The tendency of the vortex filaments to twist could be explored more in detail by examining the 
torsion, but we refrain from this calculation due to larger errors which would arise in computing third derivatives. A 
better approach would be to model the cascade more microscopically using the Gross-Pitaevskii equation, although 
the range of scales which would be available is limited. 

The possibility that the geometry of superfluid turbulence undergoes a qualitative change at low temperatures, 
where the friction can be ignored, was first considered by Tsubota and collaborators [t^ . They noticed that if a 
vortex tangle is computed in the absence of the smoothing effect of the mutual friction, its appearance is "kinky". 
Further work [t^ showed that, even in the presence of a small amount of friction, the vortex tangle is a fractal 
object. This result stimulated the study of the relation between energy and length, and between energy and fractal 
dimension [73 - l76j . The results which we present here show that, in its simplest form, the fracture nature of superfluid 
turbulence arises from the Kelvin- waves cascade along individual filaments. The actual value of the fractal dimension 
depends on the level of acoustic dissipation of kinetic energy, which in our Biot-Savart simulations is modelled by 
the numerical dissipation via the parameter S. Therefore, for decaying turbulence as in our calculations, cd depends 
on time. The fact that we observe some secondary reconnections only in the evolution shown in Fig. [50] means that 
probably fractality is a general feature of the Kelvin- waves cascade, not the result of a sequence of self-similar twisting 
and secondary reconnections; it would be interesting to pursue this aspect of the problem with more computing power 
and a smaller S. It is worth remar king that in classical fluid mechanics the relation between turbulence and fractals 
has been noticed for some time |77l - [8l| . 

Finally we notice that there are aspects of what we have discussed which relate to classical fluid mechanics, par- 
ticularly to the problem of the complexity of solutions of the Euler equation for an inviscid incompressible fiuid. In 
the absence of viscous dissipation, classical vorticity can stretch greatly, and filamentary vortex structures may twist 
wi Idly [el, [13, unable to reconnect, and may even tend to a blow-up in finite time. But there are important differences 
between incompressible Euler vortices and superfluid vortices. Superfluid vortices reconnect, lose energy by radiating 
it away, and vortex volume is not conserved (the superfluid vortex core radius is fixed by quantum mechanical con- 
straints on the rotation) . Therefore Lagrangian methods for classical flows and for superfluid flows differ in both scope 
and detail. For example, Schwarz's de-singularization of the Biot-Savart integral docs not involve a time-dependent 
core size. Another example is the cutoff S, which models sound dissipation and the reconnecting distance. In our 
case, dissipation and reconnections are essential ingredients of the problem, and the numerical method must model 
their existence. In the classical case, dissipation and reconnections are not part of the problem. 
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FIG. 1: Single vortex line initially perturbed by two Kelvin waves of wavenumbers fc = 3 and = 4 (where k = 2-k/\) and 
amplitudes A = 35 and A — 1.25, where 5 = 0.001cm. The vortex line is contained in a cubic periodic box of size D = 0.1 cm. 
This initial condition does not evolve into a Kelvin-waves cascade. 
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FIG. 2: Time evolution of three parallel vortex lines in a cubic periodic box of size D = 0.1 cm. Initially, at f = (left), 
each vortex line is perturbed by two Kelvin waves of wavenumbers fc = 3 and k — 4 and amplitudes A — 35 and A = 1.25, 
as in Fig. [1] The waves are offset by random phases. This initial condition triggers a Kelvin-waves cascade, as shown by the 
increasing wiggliness at time t — 0.0082 s (right). 
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FIG. 3: Single vortex ring of radius R = 0.1 cm perturbed by two Kelvin waves of mode numbers m = 2 and m = 5 (where 
A = 2nR/m) and amplitudes A — 5S and A — 1.45, where 5 — 0.001 cm. The vortex ring is contained in an infinite volume 
and the box (of size D — 0.1cm) is for visualization only. This initial condition does not evolve into a Kelvin-waves cascade. 




FIG. 4: Evolution of curvature and energy corresponding to the non-cascading initial condition shown in Fig. [S] Top: PDF(C) 
vs C (cm~^) at different times. Solid line: aX t = Q; dashed line: sk, t — 0.5 s; dot-dashed line: at t = 1 s. Note that the peak 
of the PDF does not move significantly to larger values of C with time t. Bottom: Energy E (cm'^/s^) vs time t (s). 
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FIG. 5: Cascading vortex configuration. Time evolution of three vortex rings of radii Ri — R2 = R3 — 0.16 cm. Each ring has 
been perturbed by two Kelvin waves. Two rings have waves with mode numbers m = 3 and m = 4 and amplitudes A — 5S and 
A = 1.45; one ring has waves with m = 2 and m = 3 and A = 5S and A — 1.45. As in Fig. [J] this initial condition triggers a 
Kelvin- waves cascade (note the increasing wiggliness of the filaments). The calculation is performed in an infinite volume: the 
box shown (for visualisation purpose) is a cube with sides of length 0.1cm. Top left: at time t = s; Top right: at t — 0.1 s; 
Bottom left: at t = 0.2 s; Bottom right: at t = 0.4 s. 
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FIG. 6: Vortex length L (cm) vs time t (s) corresponding to Fig. [S] 
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FIG. 7: Energy E {era? /s^) vs time t (s) corresponding to Fig. [S] 
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FIG. 8: Top: maximum speed Vmax (cm/s) vs time t (s) (top) and maximum acceleration {dv/dt)„-iax (cm/s^) vs t (s) (bottom) 
of the vortex points Sj {j ~ 1, - ■ ■ N) corresponding to Fig. [S] 
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FIG. 9: Curvature vs time t corresponding to Fig. [S] Top: mean curvature < C > (cm"'^) vs t (s). Middle: minimum 
curvature Cmin (cm~^) vs t (s). Bottom: maximum curvature Cmax (cm^^) vs t (s); the dashed line is the maximum curvature 
which corresponds to the numerical discretization. The maximum curvature which is observed is therefore always much less 
than the numerical resolution. 
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FIG. 10: Probability density function of the curvature PDF(C) vs C (cm~^) corresponding to Fig. [5] at selected times. Solid 
line: at i = 0. Dashed line: at i = 0.5 s. Dot-dashed line: aX t = 1 s. Note the curve initially moves greatly to the 
right (compared to the non-cascading vortices shown in Fig. |4jtop)) and then, at later times, slightly to the left (hence lower 
curvature) as the turbulence slowly decays. 
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FIG. 11: (Color online). Exponent a vs time t of the probability density function of the curvature, PDF(C) ~ C" for large C, 
corresponding to Fig. [S] The (red) dashed line is the running average. 
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FIG. 12: Analysis of the evolution of the initial vortex configuration shown in Fig. [5] (top left) in the absence of vortex 
reconnections. Top: Energy E (cm^/s^) vs time t (s). Bottom: Corresponding mean curvature < C > vs t, miniumum 
curvature Cmin vs t, and maximum curvature Cmax vs t. As in Fig. |9jbottom), the dashed line is the maximum curvature 
which corresponds to the numerical discretization. 
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FIG. 13: (Color online). Top: Smoothed vortex configuration. The (red) circles denote the actual discretization points Sj 
{j = 1, ■ ■ •A'') of the vortex configuration corresponding to Fig. [5] at time t = 0.18 s. The (blue) line is the smoothed curve 
curve Ssrnooth obtained by cubic spline fit interpolation over every n"^ — 15 points. Bottom: Slope /3 of the amplitude spectrum 
A{k) ~ for large k plotted as a function of n. Note the plateau 10 < n < 20 where the spectral slope is roughly constant. 
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FIG. 14: Mean amplitude < a > (cm) ol the Kelvin waves vs time t s corresponding to Fig. 1131 Note the saturation for 
t > 1 (s). 
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FIG. 15: (Color online) Log log plot of Kelvin wave amplitude spectrum A{k) (cm'') vs wavenumber k (cm~') at time 
t = 0.18 (s), corresponding the Fig. [TS] The (red) dashed line represents the large-k behaviour A(k) ~ fc'' with the best fit 
13 = -3.1 
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FIG. 17: (Color online). Angle between vortex filament and smoothed vortex filament. Top: Probability density function 
PDF(6) vs 9 of the angle (in degrees) between the smoothed vortex filament and the actual vortex filament ai t — 0.2 s. 
Bottom: Probability density function PDF(ze) vs zg of the scaled angle zg = {9 — 9)/a{9) where 9 and (t{9) are respectively 
the mean and the variance of the distribution of the angle 9. The dashed (red) line shows the normal PDF. Note the slight 
departure from Gaussianity at large angles. 




FIG. 18: Non-cascading reconnecting vortex rings. Time evolution of two rings of radius R = O.lGcm initially linked to each 
other. Note that the reconnection cusp relaxes and induces Kelvin waves, but there is no cascade. Top left: At time t = s; 
Top right: t = 0.16 s; Bottom left: t = 0.5 s; Bottom right: t = 0.8 s. 
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FIG. 19: Evolution of curvature and energy corresponding to Fig. 1181 Top: PDF(C) vs C (cm~^) at t = (solid line), t — 0.2 s 
(dashed line) and t = 1.6 s (dot dashed line) Note that the peak of the PDF does not move to larger values of C. Bottom: 
Corresponding behaviour of the energy E (cm^/s^) vs time t (s). 



30 




FIG. 20: Cascading reconnecting vortex rings. Time evolution of two rings of radius i? = 0.16 cm initially linked to each other. 
This initial configuration differs from the configuration of Fig. (|15|) only because some Kelvin waves have been superimposed 
to the vortex rings. Note the development of the Kelvin waves cascade. Top left: At time t = s; Top right: t = s; Bottom left: 
t = s; Bottom right: t — s. 
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FIG. 21: Top: Probability density function of the curvature PDF(C) vs C (cm^^) corresponding to Fig. (|20p at t = (solid 
line), t = 0.2 s (dashed line) and t = 0.8 s (dot dashed line) Note that the peak moves to larger values as in Fig. 1101 Bottom: 
Corresponding behaviour of the energy E (cm^/s^) vs time t (s). 
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FIG. 22: Correlation dimension cd as a function of time t (s) corresponding to Fig. 1201 Note a rapid increase of cd after the 
reconnection aX t = 0.008s, and the slow decrease which follows due to numerical dissipation. 



